Stability of the fermionic gases close to a p-wave Feshbach resonance 
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Wo study the stability of the paired fermionic p-wave suijerfiuid made out of identical atoms all 
in the same hyperfine state close to a p-wave Feshbach resonance. First we reproduce known results 
concerning the lifetime of a 3D superfluid, in particular, we show that it decays at the same rate 
as its interaction energy, thus precluding its equilibration before it decays. Then we proceed to 
study its stability in case when the superfluid is confined to 2D by means of an optical harmonic 
potential. We find that the relative stability is somewhat improved in 2D in the BCS regime, such 
that the decay rate is now slower than the appropriate interaction energy scale. The improvement in 
stability, however, is not dramatic and one probably needs to look for other mechanisms to suppress 
decay to create a long lived 2D p-wave fermionic superfluid. 

PACS numbers: 74.20.Rp, 03.75.Ss, 34.50.-s 



INTRODUCTION 



Recent success with the BEC-BCS crossover experi- 
ments in the atomic fermionic gases with s-wavc Fesh- 
bach resonances [1-3] inspired studies towards creating 
p-wave fermionic superfluids using p-wave Feshbach res- 
onances [4-7]. A number of novel features made the 
p-wavc superfluids attractive, as discussed at length in 
Ref. [8] . First of all, it is enough to put atoms into iden- 
tical hyperfine states to suppress their s-wave scatter- 
ing, leaving p-wave scattering as the strongest scattering 
channel. And indeed, p-wave Feshbach resonances be- 
tween atoms in identical hyperfine states had been iden- 
tified and studied some time ago [9, 10]. Next, the p-wave 
superfluids have a number of features distinguishing them 
from their s-wave counterparts. A richer p-wave order 
parameter allows for a possibility of observing different 
phases of the p-wavc superfluids, some of which are akin 
to the phases of superfluid Helium III [11]. Chiral and po- 
lar phases of the p-wave condensates are possible, which 
differ by the projection of the angular momentum of the 
Cooper pairs (or molecules) of the condensate onto the 
chosen axis. If that projection is ±1, the condensate is 
called chiral, and if it is 0, the condensate is called polar. 
Another important feature is that as the system is tuned 
from BCS to BEC, it does not go through a crossover as 
in the s-wave case, but rather goes through a phase tran- 
sition, as was first discussed by G. Volovik long before 
current experiments on the BCS-BEC systems became 
possible, Ref. [12]. Thus the BCS and BEC are two dis- 
tinct phases of p-wave condensates (with either of them 
possibly being chiral or polar, bringing the total num- 
ber of phases to four). Finally, when confined to 2D, the 
chiral BCS phase of the p-wave superfluids is topological 
and its vortices have trapped quasiparticles which obey 
non-Abelian statistics [13, 14]. Such quasiparticles were 
suggested to be used as topologically protected qubits to 
construct decoherence free quantum computers [15]. 



However, the program to create these superfluids suf- 
fered a setback when experimental studies of the p-wave 
Feshbach molecules showed they were unstable, with the 
lifetime varying between 2 and 20 ms [16-19]. Although 
some of these studies were done with molecules made of 
atoms of ^°K, which are inherently unstable due to dipo- 
lar relaxation [17], the rest of the studies used atoms of 
^Li, whose p-wave molecules should not, by themselves, 
exhibit any instability. 

A common mechanism which can lead to instability in 
atomic gases is the atom-molecule and molecule-molecule 
relaxation [20]. This is the process which, for example, 
involves one of the atoms approaching a molecule, with 
the result being that the molecule collapses into one of its 
strongly bound states, while the excess energy is carried 
away by the atom. Such processes are suppressed in the 
s-wave superfluid due to the Pauli principle, as was con- 
vincingly demonstrated in Rcfs. [20, 21]. However, the 
Pauli principle does not protec;t the p-wave superfluids, 
potentially leading to much shorter lifetimes. 

Refs. [22, 23] examined the stability of the p-wave con- 
densates due to these relaxation processes, as well as due 
to a possible recombination into trimcrs [22-24]. We es- 
tablished that the decay rate of a condensate of p-wave 
molecules close to Feshbach resonance is given by 
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where i is the typical interatomic spacing and is the 
van der Vaals length (the interaction range), typically 
estimated to be ~ 50 au (so the ratio I /Re, assuming 
that (. ~ 10000 au, is of the order of 200). n^/{m(,^) is 
the Fermi energy of the gas. Thus for a gas of Fermi en- 
ergy about lOKHz, this gives an estimate of the lifetime 
of l/F ~ 20ms, which is not far from what is measured 
experimentally. This should be compared with the cor- 
responding expression for the decay rate of the s-wave 
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which is orders of magnitude slower than the p-wa.ve rate, 
leading to a stable condensate. 

The calculations leading to the expression Eq. (1) were 
done solely in three dimensional space. Yet the most in- 
teresting p-wave condensate, the one with non-Abelian 
quasiparticles, has to be confined to two dimensions. The 
confinement may affect the lifetime of the condensate. In 
the absence of any experiments where the p-wave reso- 
nant gases arc confined to 2D, it is imperative that a 
theoretical calculation is done estimating this lifetime. 

In this paper we estimate the lifetime of the p-wave 
condensates close to Feshbach resonance confined to 2D. 
For purely 2D condensates we find that their decay rate 
is given by 
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This is even faster than the 3D case, Eq. (1). However, 

for the quasi- 2D condensates, the ones which are confined 
to a "pancake" of width d, where d <^ £ and at the same 
time d ^ Re, we find that the decay rate is given by 
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Notice that this expression interpolates between Eq. (1) 
and (3). Indeed, as d becomes smaller than i, Eq. (1) 
gets replaced by Eq. (4). As d is decreased, it eventually 
becomes smaller than Rg, at which point Eq. (4) gets 
replaced by Eq. (3). 

The rate in the quasi-2D geometry, given by Eq. (4) 
is somewhat faster than the 3D rate Eq. (1). So one 
may jump to the conclusion that the quasi-2D geome- 
try in fact decreases the lifetime of the condensate. This 
however must be contrasted with the fact that in 2D the 
interactions are stronger. And indeed, the typical inter- 
action energy per particle of the 3D condensate (assum- 
ing that it is in the "strong" resonance regime [22] and 
concentrating, for simplicity, on the BCS regime only) 
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that is it is of the same order as T^r, from Eq. (1). Thus 

the condensate decays in 3D as fast as it interacts. On 
the other hand, in 2D the interaction energy is given by 
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Comparing it with r2D from Eq. (3), we see that the 
interaction energy is weaker than the decay rate, by a 
logarithmic factor. 



However, the case we are interested in is quasi- 2D, 
when the condensate is confined to a pancake of width 
d. Under these conditions, the decay rate is given by 
Eq. (4), while the interaction strength is still given by 
Eq. (6), with d substituted for Rg. Provided that 
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the interaction energy in quasi-2D can be larger than the 
decay rate, thus creating a situation where the conden- 
sate might have sufficient time to form. In turn, since 
d ^ Re and logarithms, even of large arguments, are 
typically not very large, Eq. (7) may indeed hold. The 
fact that the quasi-2D BCS p-wave condensates are more 
stable than their 3D counterparts is the main conclusion 
of this paper. 

The rest of the paper is organized as follows. 

In section II we go over the analysis of the stability of 
the 3D p-wave supcrffuid, mostly following discussions in 
Ref. [22]. In particular, in subsection II A we study the 
two channel model describing the superffuid and explain 
the difference between strong and weak resonances, as 
well as narrow and wide ones, while in subsection II B we 
go over the stability analysis. 

In section III we present the analysis of the stability 
in 2D. First, in subsection III A we set up a 2D p-wave 
gas. Next, in subsection IIIB we discuss the stability of 
the condensates when the molecules are large, relevant in 
3D s-wave and 2D p-wave cases. In the next subsection 
IIIC we set up the three-body problem which needs to 
be solved to compute the decay rate. In subsection III D 
this problem is solved. Finally, in subsection HIE the 
implications of the solution are discussed and the decay 
rate is derived. Subsections III A and HID are the most 
technically involved parts of the paper and can be safely 
omitted at first reading. 

In section IV, we go over the analysis in quasi-2D, 
where the condensate is confined to a pancake geometry. 
This section is followed by Conclusions and two Appen- 
dices. 

Our final remark concerns the usage of the Planck con- 
stant h. It generally helps to omit it in calculations be- 
cause it clutters the expressions and makes them harder 
to manipulate, while it can always be restored everywhere 
by dimensional analysis. So we adopt notations where 
h = 1 everjrwhere in this paper from here on. 



II. STABILITY OF 3D p-WAVE FERMIONIC 
SUPERFLUID 



A. The two channel model 

To describe the 3D p-wave resonantly coupled super- 
fiuid we employ a two-channel model with Hamiltonian 
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Here and a arc creation and annihilation operators of 
a (spinless) fermion with mass m, while the bare spin 1 
bosonic diatomic molecule is created and annihilated by 
6^ and h^. The vector index n represents the projection 
of spin on some axis. 

The superfluid described by the Hamiltonian (8) de- 
pends on foTir parameters. Of these, co is the bare detun- 
ing, controlling the position of the Feshbach resonance, 
while the particle number N is the expectation value of 
the operator 
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Often it is convenient to trade the particle number for 
the Fermi energy tp. The Fermi energy is defined as the 
energy of a free Fermi gas whose particle number concides 
with A'' above, and is explicitly given by 
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Another way to represent the particle number is by the 
interparticle spacing 
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where n denotes the density of particles n = N/V. 

The interaction term of the Hamiltonian turns two 
fermions into a boson and vice versa. It is described by 
the interaction strength 5(|p|), the momentum depen- 
dence of which originates in the fact that the interaction 
is not point-like, rather the interaction strength is pro- 
portional to the momentum space wavefunction of the 
molecule. As such, the interaction needs to be supple- 
mented with a cut-off of order A ^ I /Re where is the 
physical size of the molecule, or the interaction range as 
discussed in the Introduction. We will write the interac- 
tion strength as 



5(|P|)=5^(|P|/A). 
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where g is the asymptotic value of the interaction 
strength for small momenta and the dimensionless func- 
tion ^ describes the quick fall-off of the interaction 
strength for momenta |p| A. In this paper, for sim- 
plicity we will let 
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with O the usual stepfunction. The precise shape of the 
cut-off function does not affect the conclusions of this 
paper, see the discussion in Section III. The parameters 
eo, cf, 5, A completely characterize the p-wave superfluid. 

Two important dimensionless combinations can be 
constructed out of these parameters. One is given by 



7' 



(14) 



If this parameter is small, 7 <C 1, then mean field the- 
ory can safely be employed to analyze the Hamiltonian 
Eq. (8). In a typical experiment 7 is indeed small, be- 
ing on the order of 7 ~ 1/10 [6, 8]. By analogy with 
s-wave Feshbach resonances, the case of small 7 can be 
termed that of narrow resonance (while the experimen- 
tally irrelevant case of 7 ^ 1 can be termed broad p-wave 
resonance). 

The second parameter. 
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can also be formed. Following Ref. [22] we term the su- 
perfluid with large C2 the case of strong p-wave resonance, 

and correspondingly the case with small C2 weak reso- 
nance. In current experiments, C2 is typically large, thus 
the resonances studied so far were strong [8]. 

It should be emphasized that under the condition C2 » 
1, it is possible to trade [8] the two-channel model (8) for 
the one-channel model 
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While C2 is indeed large in current experiments, we pre- 
fer to employ for our analysis the two-channel Eq. (8). 
Indeed, to analyze (16), one typically needs to employ 
a Hubbard-Stratonovich transformation to turn it into 
something similar to Eq. (8) first, and proceed from there. 
We find it more straightforward to work directly with 
Eq. (8). 

To further elucidate the meaning of the model Eq. (8) 
and its relation to real phenomena, we consider the scat- 
tering of two atoms with momenta k and — k into mo- 
menta k' and — k'. This scattering proceeds via the p- 
wave channel and its scattering amplitude is given by 
[28] (we emphasize that this is a partial scattering am- 
plitude, while the full amplitude is given by the standard 
expression 3/i(A:)Pi(cos 6*) where is the angle between 
the incoming and outgoing momenta) 



(17) 



Here v is called the scattering volume. In terms of the 
parameters of Eq. (8) it is given by [8] 
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In an experiment, eo (and thus loq) is tuned by varying 
the magnetic field. This induces a change in v, with l/v 
crossing zero as eg is decreased (in this regard, the be- 
havior of the scattering volume v is completely equivalent 
to the behavior of the scattering length a in an s-wave 
Feshbach resonant sc;attering). 

The parameter fco replaces the "effective range" pa- 
rameter To of the s-wave resonances and is given by 
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The poles of the scattering amplitudes describe resonant 

scattering at positive luq and bound states (molecules) at 
negative wq- These occur at (neglecting the ik^ term in 
the denominator of Eq. (17), small at small ojq) 



fc2 
m 



This elucidates the meaning of wq introduced in Eq. (18). 
One can also remark that the total scattering cross sec- 
tion, for the p-wave scattering, is given by [28] 



a = 127r|/i(fc)r 



(21) 



ujQ, as long as it is positive, can be measured in an ex- 
periment by looking at the energy of colliding particles 
at which the scattering cross section a has a maximum. 
At the same time, eo is varied by the magnetic field ac- 
cording to 
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where Bo is the magnetic field corresponding to the res- 
onance and fj,B is the effective Bohr magneton. Thus 
measuring lvq and B — Bo simultaneously allows to de- 
termine whether the resonance is weak or strong (c2 <C 1 
or C2 ^ 1). Strong resonances will appear as the ones 
where the slope of the curve lvq vs hb {B — Bq) is small. 

The dependence of the scattering volume v on the mag- 
netic field B is given by 
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reminiscent of the magnetic field dependence of the scat- 
tering length in an s-wave Feshbach resonance experi- 
ment. We notice that if one tunes the magnetic field 
B off resonance by the amount such that ^i{B — Bq) is 
equal to the Fermi energy of the gas, then the scattering 
volume will be much smaller than the cube of the inter- 
particle spacing if the resonance is narrow (7 <C 1), 
and much larger than the cube of the spacing if the reso- 
nance is broad (7 ^ 1). This justifies the name "narrow" 
vs "broad" . Indeed, in case of the broad resonance, v vs 
B — Bo graph will appear much broader than in the case 
of narrow resonance. Notice the complete equivalence of 
this term with its s-wave counterpart [6, 8, 29-33]. 



We emphasize that other authors, such as the ones of 
Ref. [23] , prefer to restrict the usage of the term "broad" 
vs "narrow" for C2 being large or small. While the con- 
crete terminology is a matter of taste, 7 is a parameter 
which more accurately reflects the notion of "broad" vs 
"narrow" , as these terms are used in the s-wave resonance 
context. 



B. The analysis of the stability of the 3D p-wave 
condensate 

Here we reproduce the analysis of the stability of the 
3D p-wave condensate from Ref. [22]. Suppose, at some 
negative luq. p-wave molecules form whose binding energy 
is LUo- The radial part of their wavefunction, at distances 
much larger than the interaction range Rg ^ 1/A, is given 



(20) by [18, 28] 
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where n = \Jra jciJol- At distances much smaller than 
1/k but still much larger than R^, Re <^ r <^ 1/ k, the 
wavefunction can be well approximated by 
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We notice that the normalization condition of the wave- 
function 
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leads to the following normalized expression 
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(27) 



independent of luo- In other words, most of the weight of 
the wavefunction is concentrated at distances Re, or the 
molecules are small. 

Suppose two such molecules collide. It is possible that 
the collision will lead to one molecule forming a strongly 
bound state, while the atoms of the other molecules ab- 
sorb the energy and fly apart. It is also possible that 
three of the atoms form a strongly bound trimer, while 
the remaining atom flies away [22, 23] . The rate of this 
process can be estimated as follows. The total rate is 
given by 
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where n = 1/^^ is the density of particles, v is their veloc- 
ity, and (7 is the inelastic cross section for this process. In 
turn, a can be estimated as a product of the elastic cross 
section of two small objects size Rg each, or R^ times 
the time they spend in the vicinity Re/v times the rate 
of the decay. That rate may be hard to calculate, but it 
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must be of the order \/{mR^) by dimensional analysis. 
Putting this all together we find 
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which matches Eq. (1) given in the Introduction. 

Another simple way to derive F proceeds as follows. 
Due to the ultraviolet divergencies in the p-wave two- 
channel model Eq. (8), any transition amplitude, elastic 
or inelastic, between some states of two molecules or a 
molecule and an atom, should go as (up to a phase). 
The inelastic cross section can then be estimated as [281 
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where kf = I /Re is the final momentum of the particle, 
and ki = mv where v is the incident velocity. This leads 
to 
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which coincides with Eq. (29). This derivation makes it 
clear that Eq. (29) applies not only to molecule-molecule, 
but also molecule-atom collisions. 

On the BCS side of the resonance, one could argue 

that the atoms spend some fraction of their time virtu- 
ally forming molecules. Those will also undergo atom- 
molecule relaxation so the decay rate Eq. (29) applies 
here as well. One can also consider direct 3 body recom- 
bination of atoms. Experimcnital and theoretical studies 
[9, 34] lead to a rate numeric;ally close to the one esti- 
mated here for the molecule-molecule and molecule-atom 
relaxation, at the relevant densities. 

Is the decay rate given by Eq. (29) too fast, or suf- 
ficiently slow? Obviously it has to be slow enough for 
the condensate to form and equilibrate before it decays. 
Issues of equilibration are subtle. Instead of consider- 
ing them, we estimate the interaction rate between the 
atoms in the BCS phase and between the molecules in the 
BEG phase. This rate corresponds to a time scale which 
is clearly shorter than the equilibration time. Thus, if 
this energy scale is larger than the decay rate, there is 
a chance that the condensate would have time to equili- 
brate. 

In the BCS regime, the molecules have positive energy 
and if left in vacuum decay into atoms. The molecular 
decay rate, computed at the Fermi energy, corresponds 
to the time scale at which atoms interact. This rate can 
be read off the scattering ampHtude Eq. (17). For that, 
we find the pole of this scattering amplitude, not just the 
real part as in Eq. (20), but also with its imaginary part. 
It is given by 
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We estimate it in case of strong resonances, or C2 S> 1. 
Substituting C2 from its definition Eq. (15) and choosing 
Wo to be the Fermi energy of the gas, or Wq ~ l/{mi^), 
we find the interaction energy 



1 Re 
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Note that the ratio g'^ /{I + C2) is an increasing func- 
tion of C2. Thus for weak resonances when C2 -C 1, the 
interaction energy would be even weaker than the one 
given by Eq. (33). So the case of the strong resonance 
corresponds to the strongest possible interactions. This, 
combined with the fact that the p-wave Feshbach reso- 
nances experimentally studied so far appear to be strong, 
prompts us to concentrate in this paper mostly on strong 
resonance. 

Yet as we see this interaction has the same functional 
form as the decay rate Eq. (29). Of course, these are 
estimates of these quantities and perhaps numerical co- 
efficients omitted in Eqs. (33) and (29) conspire to make 
one larger than the other. But typically we expect that 
the decay rate and the interaction energy are of the same 
order, precluding the formation of the condensate before 
it decays even in case of strong resonance. 

Another way to estimate the interaction rate is in the 
BEC regime. With the elastic scattering cross section 
being R^, the elastic scattering rate is 
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since mv is the characteristic momentum of the 
molecules, which is roughly equal to l/£. Obviously, 
Re/i "C 1, thus the interaction rate between the 
molecules in the BEC regime is even slower than the 
interaction rate between the atoms in the BCS regime. 
In particular, it is slower than the decay rate of the 
molecules, thus making the observation of the 3D BEC 
p-wave condensate even less likely than its BCS counter- 
part. 

We can compare the estimates derived here with the 
measured decay constants from Ref. [19] (see their Table 
I). The decay constant is defined as 



(35) 



and for the estimate Eq. (29) gives (for once, we explicitly 
reintroduce the constant h) 



K = 
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(36) 



The imaginary part is the decay rate. 



Here we take Re = 50 au. and take m to be the mass of 
^Li. This is very close to the measured atom-molecule de- 
cay constant and is one order of magnitude smaller than 
the measured molecule-molecule decay constant. We do 
not know why the measured molecule-molecule decay 
constant is faster by a factor of 10, but note that the 
derivation presented here ignores the details of the short 
range physics and could easily be off by a factor of 10. 
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Wc also note that Rcf. [19] quotes that the elastic scat- 
tering rate between the molecules is faster than the in- 
elastic rate, while the estimates presented here point to- 
wards the elastic rate being slower than the inelastic rate. 
We do not know the reason for this discrepancy. 

The conclusion is, the 3D p-wave superfluids decay as 
fast as they interact and do not have time to form before 
they decay. 



fermionic loop separating molecules of spin n and p is 
diagonal in spin indices and is denoted Yi^i, = The 
full propagator is then given by 



DO{p) - S(p) 



D{p)6^^, 



(39) 



with the fermionic loop taking the value (a factor 2 ap- 
pears from indistinguishability of the fermions) 



III. THE p-WAVE SUPERFLUID IN TWO 
DIMENSIONS 

A. The two channel model 

We now consider the case of the 2D j5-wave superfluid, 
governed by the same two-channel model Eq. (8), but in 
two dimensional space. Reducing the space dimensional- 
ity changes the dimension of the coupling. In fact, the 
p-wave two channel model's upper critical dimension is 2 
[35, 36]. The coupling constant g is now a dimensionless 
quantity. The linear divergence which led to the appear- 
ance of C2 in the 3D calculations is now replaced by a 
logarithmic divergence. 

As in 3D, it is instructive to compute the scattering 
amplitude of two atoms. Since to our knowledge this 
was not done before in the literature for the two-channel 
model, here is its derivation. 




FIG. 1: The rcnormalizod propagator of spin 1 molecules. 
Here the standard notations for the two-channel model Eq. (8) 
are used: the thin straight lines are fermionic propagators, the 
thin wavy lines are bare bosonic propagators, and the thick 
wavy lines are full bosonic propagators. 

According to the Hamiltonian (8) the propagator of 
fermionic atoms is the free propagator 
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Po — P^/'itn + iO' 
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For simplicity of notation, p is used both as the three- 
vector (p,j'o) and as the absolute value of the momentum 
|p|. The bare propagator of the bosonic spin 1 molecules 
is equal to 



Po -pP'llm- eo + iO 



- = D\p)b^,. (38) 



The molecular propagator should be renormalized by the 
presence of fermionic loops as depicted in Fig. 1. The 
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All singularities in the complex pq plane lie slightly below 
the real axis. The molecular propagator becomes 

D{p) = 



Po 



I^^-£o + c(Po - 1^) log 



1 + 



Here, 



eo 



m.g^A^ 
47r ■ 



is a renormalized detuning, while 
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(41) 



(42) 



(43) 



is the constant controlling the strength of the Feshbach 
resonance, equivalent to C2 in the 3D calculations. 

Just as in 3D, the molecular propagator has a pole at 
p = and when po is taJsen to an appropriate value. 
We denote the real part of the value of pq at the pole 
as ojQ. The pole corresponds to the binding energy of 
the molecule if wo < (in which case the pole occurs 
at Pq = Wo) and to the resonance if > (in which 
case only the real part of po is equal to at the pole, 
while the imaginary part of po describes the decay rate 
of positive binding energy molecules in free space). The 
value of Wo is controlled by tuning eo. 

We can calculate wq from the condition that it is the 
real part of po at the pole of the molecular propagator. 
This gives 



e'o = wq + cwolog 



A2 



mu)o 



(44) 



Then the (physical) molecular propagator is 



7 



D{p,po + uJo) = 
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Notice that for future convenience we shift the energy po 
by Wo so that po in Eq. (45) measures energy from the 
binding energy. 




FIG. 2: The diagram corresponding to the scattering of two 
atoms. 

Now we can compute the scattering amplitude of two 
atoms in 2D, just like we did it in 3D in Eq. (17). The 
scattering of two atoms with incoming momenta k and 
— k into momenta k' and — k' proceeds via formation of 
a molecule, see Fig. 2. Thus the scattering two-atom T- 
matrix coincides with the propagator Dfj^^ computed at 
momentum p = 0, and at energy pg + wq = jm^ con- 
tracted with the incoming momentum kp, and the outgo- 
ing momentum /c^. The relationship between the scat- 
tering T-matrix and the scattering amplitude depends 
on the dimensionality of space. In 2D it is given by (see 
Appendix A for a derivation) 



/ = -: 



(46) 



Putting this all together gives for the partial amplitude 
of p-wave scattering, defined in Eq. (B2), 



m\fk ( eg 



m 



clog 



A2 



(47) 

Here eg can be substituted in terms of cjg using Eq. (44), 
and it is assumed that A^/fc^ ^ 1. This expression con- 
forms to the general form of the 2D scattering amplitude 
Eq. (B4) . We emphasize that even though it was derived 
from the 2D version of the two-channel model Eq. (8), 
only the parameters eg and c follow from that model. 
Other than that, any p-wave two dimensional scattering 
at low energy must take this form, regardless of the model 
used (a similar point for the 3D scattering was empha- 
sized in Ref. [8]). This argument will become important 
in section IV. 

The scattering amplitude Eq. (47) has a pole at 
fc^/m = wg at negative wg and at real part of fc^/m equal 

to loq at positive ljq, just as its 3D counterpart and as 
follows from the properties of the molecular propagator. 



B. Stability of condensates with large molecules 

The radial part of the wave function of a 2D p-wave 
molecule with energy close to zero is given by 



*(r) 



(48) 



To see this, compare with the 3D case described by 
Eqs. (24) and (25) and recall that the zero energy so- 
lutions of the Schrodinger equation go as l/r'"*"^ in 3D 
and 1/r' in 2D, where I is the angular momentum. The 
normalization condition now follows from 



rdr\"^{r)\ = log 



kRp 



(49) 



where k = ^/rn\uJo\, and where, as before, wg < is 
the binding energy of the molecule. This integral is now 
divergent logarithmically at both lower and upper limits, 
so the molecular weight is equally distributed between 
and 1/k. So the size of the molecule is no longer Rg, like 
it was in 2D, but rather 1/k. 

To compute the rate of the atom-molecule relaxation, 
we go through the same steps as we did in section II B. 
However, we need to take into account that the molecules 
no longer have size i?e, but rather 1/k. Let us do it, for 
generality reasons, in an arbitrary number of dimensions 
d. 

The decay rate is still given by 



r ^ nav, 



(50) 



where n is the density, a is the inelastic cross section, 
and V is the velocity, just as in 3D, Eq. (28). However, 
the meaning of the inelastic scattering cross section a is 
now different. We now have two objects of the size 1/k 
colliding inelastically. The inelastic cross section is given 
by the product of their elastic cross section, proportional 
to 1/k'^~^, the time the molecules spent together, given 
by 1/{kv), the collapse rate during that time, l/{mRl). 
This should still be multiplied by the probability that 
the three atoms out of four which constitute 2 molecules 
find themselves at distance i?e from each other, so that 
they were all at distances of the order of the force range 
between them. To find this probability is not a simple 
problem. We have three fermions, each pair interacting 
strongly since they are close to a p-wave Feshbach reso- 
nance. 

Suppose their 3-body wavefunction is given by the fol- 
lowing scaling ansatz 



(51) 
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Here r denotes a collective coordinate of the three par- 
ticles, such as, for example, a hyperspherical radius. 
Eq. (51) serves as a definition of a scaling exponent 7. 
Then the probability that three fermions find themselves 
at a distance Re from each other is given by 



P - (i?eK 



(52) 



The power of 2d comes about because of the phase vol- 
ume of putting two fermions at a distance Re from the 
third one, while 27 arises from the behavior of the square 



of the wavefunction, \ip\ 
together gives 



r27 



Putting all these factors 



r ~ n 



d-1 



1 1 



{ReK) 



mRl KV 

{ReKf'^^'-^ 



(53) 



Finally, close to Feshbach resonance, the size of the 
molecules is close to their separation, or 1/ac Then 
the Eq. (53) simplifies to give 



1 



2d+27-2 



(54) 



This is the final answer for the decay rate of large 
molecules. 

Let us check that this expression indeed gives the cor- 
rect answer in the case of 3D s-wave molecules (which 
are large). In this case, d = 3, and 7 « —0.22 [20, 21]. 
Then 



1 

^2 



3.55 



(55) 



as was indeed derived in Ref. [20], and as was discussed 
in the Introduction Eq. (2). For easier comparison with 
Ref. [20], recall that it is sometimes beneficial to intro- 
duce the decay constant K = V /n introduced in Eq. (35). 
Recall also that in the 3D s-wave problem, k = 1/a, 
where a is the scattering length. This gives 



K 



Re 

m 



2.55 



(56) 



the form discussed in Ref. [20]. 

We also briefly examine the case of 3D s-wave bosons 
close to Feshbach resonance. Then 7 = —2 (as a con- 
sequence of the presence of Efimov states) and Eq. (54) 



gives 



1 



(57) 



A decay constant for the boson problem is defined as 
K = r/nP. Substituting the scattering length a for i we 
find 



K 



a 

m ' 



(58) 



a well known form for the boson problem [37-40]. 

When apphed to the 2D p-wave problem, Eq. (54) gives 



1 



Re 



2-1-27 



(59) 



At issue now is calculating the exponent 7. This calcu- 
lation is the subject of the next two subsections. 



C. The three-body problem 

We need to compute the three body wavefunction close 
to Feshbach resonance. This can be done in either co- 
ordinate or momentum space. We are going to present 
the momentum space derivation, since it can be done 
using the standard techniques of many-body theory. 




FIG. 3: The diagram corresponding to the wave function of 
free fermions. The square block represents the atom-molocule 
T-matrix, and the three outgoing fermionic lines represent the 
three fermions whose wave function we would like to compute. 

First of all, let us show that the scaling of the three 
body wave function is related to the behavior of the 3- 
body scattering amplitude. Suppose a three body scat- 
tering matrix T, depicted in Fig. 3, is known. It is a 
function of the incoming momenta and energy, and the 
outgoing momenta and energy. To arrive at the wave 
function of three particles, we fix the incoming momenta 
and energy to be on shell. Then we multiply the T ma- 
trix by the four outgoing Green's functions, one bosonic 
and three fermionic. Finally we Fourier transform with 
respect to the fermionic outgoing momenta, respecting 
the momentum-energy conservation. For the illustration 
of this procedure, see Eq. (Al) which represents this in 
case of just one particle scattering off a potential. 

Suppose the T-matrix scales as , where p represents 
the overall momentum scale. Then it is easy to figure out 
the scaling of the wave function. There are three outgoing 
energies and momenta, corresponding to the three out- 
going fermionic lines. However, the energy-momentum 
conservation restricts the number of linear independent 
energies and momenta to two. Hence there are two inte- 
grals over energy and momentum, contributing the power 
2(d-|-2) (energy is counted as momentum squared). There 
are four outgoing propagators, one bosonic Eq. (45) and 
three fermionic Eq. (37), contributing the power —8 (the 
bosonic propagator includes logarithms, but these are ir- 
relevant for the purpose of power counting) . Finally there 
is a vertex where a bosonic line splits into two fermions, 
contributing a power 1 due to the p-wave momentum 
dependent factor. Putting it all together, and remem- 
bering that the coordinate scaling is opposite in sign to 
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the momentum scaling, we find that the scaUng of the 
wave function is given by 

7=-(7' + 2((i+2)-8+l), (60) 

or in two dimensions, 

7 = -y - 1- (61) 

To proceed, we need to know the scaling of the T- 
matrix, representing the scattering of a fermion and a 
spin 1 molecule. 

D. Solution of the three body problem 

In the scattering problem, let the incoming molecule 
have spin /i and the outgoing spin The T-matrix will 
then in general be a tensor T^,,^. In the center of mass 
frame, with the incoming molecule having momentum 
k and the outgoing p, the tensor T^^, consists of five 
terms proportional to p^Pi,, p^fc;y, k^p,^, and fc^fc^. 
However, we are interested in the wave function at short 
range r <C which corresponds to large p, or p ^ k ^ 
K. Then it is sufficient to set k = 0. Thus define [24] 

7).i.(p,Po) = Ti{p,po)6^^ + T2{p,po)p^Pu/p'^ 

= ^ r,(p,poX(p). (62) 

i=l,2 




FIG. 4: The series of diagrams which sum to give the scatter- 
ing T-matrix of a spin 1 molecule and a fermionic atom. 

The scattering T-matrix is the sum of the series of di- 
agrams indicated in Fig. 4. These diagrams appear the 
same as those studied in the s-wave three-body problem 
[41, 42], however the Feynman rules are of course dif- 
ferent in the present problem. For weak resonances, the 
diagrams form a perturbative series in which only the 
first few diagrams may be kept. However, for the strong 
resonances studied in the present paper the diagrams are 



all of the same order and the sum of the series of diagrams 
is needed. The summation may be attained by construct- 
ing a Lippmann-Schwinger type integral equation for the 
scattering amplitude. This summation was performed in 
a similar manner in the three-dimensional problem stud- 
ied in Ref. [22] (see also Ref. [23]). The kincmatic;s 
are chosen as follows; the incoming molecule is on-shell 
with three-momentum (0,a;o) while the incoming atom 
has three-momentum (0,0). The outgoing molecule has 
(P)Po + Wo) and the outgoing fermion (— p, — Po)- The 
scattering matrix does not include external lines. Then 
the integral equation takes the form 

T^,^{P,P^ = g{p)g{p/2)G{p,po + u;o)p^,Pu 
-4« j "(^;^^Ma(q> qo)G{q, -qo)G{p + q,po + 90 + wq) 

xZ)(g,go + t^o)(p+ |)a(g+ |).5(|P+ |l)5(|q+ |l)- 

(63) 

Repeated indices are summed over. Z is the residue of the 
molecular propagator at the pole and is needed for correct 

normalization of the scattering matrix. It is a function 
of Wo and c whose precise value will not be needed in the 
following. 

The integral over qq in Eq. (63) may be performed 
by closing the contour in the upper half plane, setting 
Qo —q^/lm. In order to solve the integral equation 
it is then convenient to let po —p^ /2m. This ensures 
that the frequency dependence of T^^ is the same on both 
sides of the integral equation. The integral equation is 
then solved for Tj(p) = Tj(p, — p^/2m). Subsequently, 
this solution can then be used to find Ti{p,po) at any 
Po < 0. 

To project onto the functions Ti and T2 defined in Eq. 
(62) multiply the integral equation (63) by 'U^^(p). The 
left hand side will then contain the matrix 

Uki = <(p)«i..(p) =(l\) ■ (64) 

V / ki 

This matrix is invertible and it is thus possible to find a 
set of coupled integral equations for the functions Ti (p) 
and T2{p). Upon multiplying Eq. (63) by U'^^ul^ip) it 
becomes 



Tj{p) = -2Zg{\p\)g{\p\/2)G{p,-py2m + u;o)p^52j j qdqD{q,-qy2m + iVo)bji{p,q)T,{q,-qy2m). (65) 

The dimensionless matrix b is given by 

with 9 the angle between p and q. 
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To further simplify, assume that both p and q are cut off at A using the usual step function. Then the matrix b 
may be calculated with the result 



611 

b2i 

b22 



\J(uj —jp' — — p^q^ +UJ —p^ — q^ 



-TT + TT 



(2a;-2p2_g2)(^_p2_g2) 



p2q2 



+ 



^(y) ~ p^ — q^)'^ — p'^q^ 
2 2 ^(uj^p 



{2uj-2p^ -q^){uj-p^ -q^) 



2\2 



p2q2 



- 2a; + 2/ + q^ 



27r 



n_ (4a;-3p^-2g^)(a;-p2-g^) 
■ 2 p^q^ 



\/ {to — pi^ — q^Y ~ V^l^ + to — p^ — q^ 



+ - 



7r/2 



^/{uJ — p"^ — q^y — p^cp' 
Here a dimensionless detuning is defined as a; = wq 



-2^^ hA. ^ + 3a; - 3^^ - 



p2q2 



(67) 

(68) 
(69) 

(70) 



According to Eq. (51), in order to determine the 

behavior of the three-body wavefunction, the functions 
Ti{p) and T2{p) are needed for momenta ^—mujQ <C p <C 
A. We now proceed to solve Eq. (65) analytically in this 
range of momenta. In the following, the limit c ^ 00 will 
be taken to ensure that the Feshbach resonances studied 
are strong. For momenta in the range of interest the 
integral equation reduces to 

dq bji{p/q), 



Tj{p)=2Zg''m52j 



^ log (f ) 
The matrix b takes the limiting forms 
K<p<g< A: 



Tiiq). (71) 



6 = 



TT I 37r p 
■ 2 8 92 
_97r 
\ 16 

K<g<Cp< A: 





Eq. (71) suggests a solution of the form 

'A^ 



T,{p) = Cilog' 



P 



T2{p)^C2l0g" 



(72) 



(73) 



(74) 



Substituting into Eq. (71) and keeping leading terms 
results in the set of linear equations for the coefficients 
Ci and C2 

Cilog"f^ 
VP 



Calog" 



+C2 



4 

3a 
4 

"3a 

log" 



Cilog" ( ^ ) +C2l0g" 



2Ci 

A 

K 



log 

- log" 



log° 



p. 

2Zg'^m. 
(75) 



This set of equations has solutions only if a = 
Matching coefficients finally results in the solutions 



[4 /A 
Ti{p) = pcos -log log ( - 

r A 

T2{p) = \/2pcos 
with the amplitude satisfying 



37r 
T 



"4 




-log log 









(76) 



(77) 



cos [| log log (A/k)] 

This solution contains a free parameter, as it does not 
allow the determination of both p and cj) independently. 

In solving Eq. (65) numerically, it is found that the 
overall amplitude of the solutions Ti{p) and T2{p) con- 
verges very slowly. We attribute this slow convergence 
to the logarithmic behavior of the solutions; in the nu- 
merical study it is important to keep log(A/K) 3> 1 while 
simultaneously the configuration space must contain a 
large number of momenta p for which log(K/p) ^ 1. Hav- 
ing this in mind, we therefore write the solutions as 



Ti{p) = Pi cos 



"4 




-log log 









= 1,2 



(78) 

rather than the separate 



and determine the ratio P2/ pi 
values of these amplitudes. 

The solutions Tiip) and 12 (p) of the coupled integral 
equations (65) are shown in Fig. 5 for a large value of 
C2 . Also shown are the analytical solutions (78) with the 
parameters 

</>! = 2.80 4>2 = .56 
pi = .71 p2 = .96 . (79) 

We observe that ^1 — (^^2 ~ 37r/4 and P2/P1 ~ V2, both 
with a 5% error. Thus we conclude that the analytical 
expressions given in Eq. (76) are correct. 
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in the 3D case, we can term the case when this param- 
eter is large as strong resonance, although unlike in 3D, 
any resonance when the gas is sufficiently dilute becomes 
strong. If the resonance is weak, then c <C 1, since the 
logarithm in Eq. (82) is always large. 

Assuming that the resonance is strong, we find the pole 
of Eq. (47) to be at 



fc2 
m 



■KUJQ 



log 



A2 



(83) 



In the opposite case of weak resonance, we find 

~ 0^0 — ITTCLOq. 



(84) 



Substituting loq = ep, we estimate from here the atomic 
interaction energy as 



FIG. 5: Ti{p) and T2{p) obtained by solving Eq. (65) for 
ujQ = —10^''^ and C2 = 10® (blue, thick). The solution uses 
Gaussian-Legendre quadrature for the numerical integration 
[43] with 2000 grid points. Also shown are the analytical 
solutions (78) with parameters chosen as in Eq. (79) (black, 
dashed) . 



E. The lifetime and the interaction energy of the 
2D p-wave condensates 

We are now in a position to finish the calculation of 
the lifetime of the 2D p-wave condensate. The scattering 
amplitude is given by Eq. (62), where in turn Ti and T2 
are given by Eq. (76). These expressions do not scale at 
all, corresponding to 7' = 0. Thus, according to Eq. (61), 



7 = -7' - 1 = -1. 
Substituting this into Eq. (59) we find 



(80) 



(81) 



which gives the decay rate of the 2D p-wave superfluid 
as discussed in the Introduction, in Eq. (3) . 

This should be compared with the interaction rate for 
the 2D p-wave atoms at the Fermi energy of the gas. 
The find this energy first the BCS regime by examining 
the scattering amplitude of two atoms at positive as 
found earlier, Eq. (47). The imaginary part of the pole 
of this scattering amplitude, computed at cjo = ef, gives 
us the needed interaction rate, quite analogously with 
Eq. (32) which we employed in 3D. The answer crucially 
depends on whether the parameter 



clog 



R„ 



(82) 



is large or small (as before, I is the interparticle separa- 
tion and Rf, = 1/A is the interaction length scale). As 



E. 



1 



2D 



log 



(85) 



in the case when the resonance is strong and 



E2D 



c 



(86) 



when the resonance is weak. We expect that in a typ- 
ical experiment the gas will be dilute so the resonance 
will be strong, and that's why we quoted Eq. (85) in the 
Introduction, Eq. (6). In either of these two cases, we 
see that the interaction energy is much smaller than the 
decay rate, given by Eq. (81), so we expect purely 2D 
p-wave gases to be unstable. 

Likewise, in the BEC regime, we need to estimate the 
elastic scattering rate of two molecules. The molecules 
(unlike atoms) scatter in the s-wave channel. The scat- 
tering of two molecules in 2D must proceed according 
to the standard rules of quantum mechanics, which pre- 
dicts that the s-wave scattering cross section of particles 
at sufficiently low momentum k goes as [28] 



1 



k log2 ( 



k 



(87) 



where the unknown constant a depends on the details of 
the interactions. This gives for the elastic collision rate 
(taking k of the order of 



E2D,BEC = O-einV ■ 



1 



mP log 



V k 



This energy is somewhat smaller than the BCS estimate 
Eq. (85), but since the logarithmic factor is not very 
large, we can think of it as being of the same order as 
Eq. (85). Thus the BEC 2D superfluid is as unstable as 
its BCS counterpart. 
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IV. QUASI-2D SUPERFLUID 

Now let us consider the last remaining question, the 
stability of a 2D p-wavc supcrfluid confined to a pancake 
of width d. Such a superfluid is still described by Eq. (8), 
but with a presence of an extra confining potential in the 
third direction, which we denote z 



V{z) = -mcj'^z^, 



(89) 



where iv is the confining frequency. The width of the 
pancake d is related to the oscillator frequency via 



d' 



The quasi-2D regime is 



(90) 



(91) 



Re ^ d because otherwise the physics of the Feshbach 
resonance is modified by the confinement (also, Rg is very 
short so that the confinement to the scale below i?e is not 
currently technologically possible), while d <C in order 
for the gas to be truly confined to 2D. 

The scattering of identical fermions close to a p-wave 
Feshbach resonance follows from this Hamiltonian. Cal- 
culating it in the confined geometry is an involved prob- 
lem, first solved for the case of s-wave resonance confined 
to 2D in Ref. [44, 45] (see also the first calculation of this 
type, done for the s-wave gas confined to ID, in Ref. [46]). 
For the case of p-wave resonances, this problem was stud- 
ied in Ref. [47] in both ID and 2D. Yet, for our purposes, 
we do not need to know the answer. It is enough to know 
that the scattering is still described by Eq. (47), albeit 
with coefficients c and Cq no longer related to the parame- 
ters of the Hamiltonian Eq. (8) the way they were before, 
but rather being some more complicated functions, which 
also depend on d among other parameters. This is be- 
cause any low energy p-wave scattering in two dimensions 
must be described by the expression Eq. (B4), with the 
function .gi(fc) having the low k expansion as in Eq. (47). 

Thus the interaction energy of the two atoms confined 
to this geometry, derived in the previous section by using 
Eq. (47) only, is still given by Eq. (85) or Eq. (86). The 
only difference is that Re in these relations should be 
traded for rf, as d is now the smallest lengthscale at which 
the 2D physics is still at work. 

Yet the interaction rate of the molecules in the BEC 
regime will actually be given by the 3D formula Eq. (34). 
This is related to the fact that in quasi- 2D geometry, the 
coefficient a in Eq. (88) is typically anomalously large, 
and leads to [45] 



E. 



quasi-2D,BEC 



^R^ 



(92) 



This energy scale is very small, thus there is no hope to 
observe the BEC of p-wave molecules, even in the quasi- 
2D geometry. So we concentrate on the case of the BCS 
phase. 



Now the decay rate in quasi-2D can be deduced in the 
following way. Atoms decay when three of them approach 
each other at distance R^. This distance is much shorter 
than d, so that the 3D decay physics must take over. So 
we should not use Eq. (81) to compute the decay rate. 
Rather, we need to revert back to the appropriate expres- 
sion in 3D, given by Eq. (29), with one modification. In 
Eq. (29), £ denotes average distance between the parti- 
cles in three dimensions, while Eq. (85) and Eq. (86) are 
written in terms of average distance between particles in 
two dimensions. These are related by the obvious 



This leads to the decay rate 



- quasi— 2D 



1 Re 

mP d ' 



(93) 



(94) 



where (. is now the two dimensional distance. This con- 
cludes the derivation of Eq. (4). 

We therefore see that the necessary condition for the 
existence of the stable superfluid in the quasi-2D geome- 
try is given by E2D > Tquasi-2D or 



max I log 



d 



l\ d 



(95) 



Since logarithms, even of large arguments, are typic;ally 
not very large (we assume c is a constant generally of the 
order of 1, although its value is controlled purely by a 
particular Feshbach resonance and it can be calculated 
from its physics on a case by case basis), it is possible 
that this condition will be satisfied in experiments. 

Indeed, a typical value for d would be, perhaps, half a 
wavelength of the light used to create a confining poten- 
tial. This gives d ~ 250 nm, or 5 • 10^ au. With £ ~ 10^ 
au, the ratio l/d ^ 2, and its logarithm is basically 1. At 
the same time d/Re can be kept just marginally smaller 
than 200. So Eq. (95) is satisfied with a large margin. 

Yet we must remember than Eq. (95) is but a neces- 
sary condition for the decay being slow enough to allow 
for eriuilibration of the superfluid. In prac;tic;e, the equili- 
bration involves many collisions between atoms and may 
go much slower than E2d- In truth we have only one 
example of the decay slow enough that that we know 
from experiment that there is enough time for equilibra- 
tion and meaningful experiments in the superfluid phase, 
that of the 3D s-wave case. In that case, the ratio of the 
Fermi energy (as a crude estimate of the characteristic 
energy of the condensate) to the decay rate is an impres- 
sive ie/Re) - 10^ The p-wave superfluids conflned 
to 2D are far from being that stable. 
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V. CONCLUSIONS 

This concludes our studies of the stability of the 

ferniionic paired supcrfluids close to p-wave Feshbach res- 
onance. From the analysis of this paper, it is clear that 
the 3D p-wave gases are inherently unstable. The situa- 
tion improves if they are confined to 2D. Yet it is nowhere 
near the case of 3D s-wave superfluids, which are sta- 
ble for all practical purposes due to a very high scaling 
power in their decay rate Eq. (55). A promising route 
to increase stability further in a quasi-2D setting seems 
to be a creative application of the optical lattice, simi- 
lar to what was recently done in Ref. [48] for a different 
problem. This should be a subject of further work. 



APPENDIX A: RELATIONSHIP BETWEEN THE 
T-MATRIX AND THE SCATTERING 
AMPLITUDE IN 2D 

Consider a particle of mass m scattering in a potential 
in 2D. Its wave function is given by 



V'(r) = e"'" + 



( ll 
J (27r: 



P_ 



E- 



2m 



T(k,p) e'P^ (Al) 



where T(k, p) is the scattering T-matrix computed be- 
tween momenta k and p at energy E = k"^ /{2m). We 
compare this expression with the definition of the scat- 
tering amplitude in 2D, given by the large r expression 
of the wave function [28] 
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\/-i = exp (-i7r/4) , (A2) 



where if is the angle between the incoming momentum k 

and the position vector r. Doing the angular part of the 
integral in Eq. (Al) at large r by the steepest descend 
method we find 



V'(r) = e*''"' + 




+ T(k,-p,) 



2tt 



-ipr 



2m 



2m 



iO 



(A3) 



where p^ denotes a vector whose length is p, but which is 
directed along r. Change variables in the second integral 
to get 



APPENDIX B: CONSTRAINTS PLACED ON 
THE SCATTERING AMPLITUDE IN 2D BY 
UNITARITY 



^"^T(k,p.),/^e--^:^A , 



47r2 



(A4) 



where the contour of integration goes above the p = 
singularity. Doing the integral by residues gives 



\Phdkr 



r(k,k^)e' 



(A5) 



Comparing with Eq. (A2) gives 



/=- 



m 



:r(k,k'). 



(A6) 



When comparing this expression with the one used in 
the text, Eq. (46), one needs to remember that m in 
Eq. (A6) is the reduced mass of two fermions and should 
be replaced according to to ^ to/2. 



The scattering amplitude in 2D is constrained by uni- 
tary, just like its 3D counterpart. Here we reproduce the 
appropriate derivation from Ref. [28]. 

Consider the scattering of a particle of mass m with 
momentum k into the momentum k' such that k = k' , 
but the angle between these two vectors is (f. The scat- 
tering amplitude can be expressed in terms of the phase 
shifts Si according to [28] 



l—oo 

iv27rK , 



Introduce the partial scattering amplitudes 



Since 



le^^^f = 1, 



(Bl) 



(B2) 



(B3) 
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a constraint on the form /; (k) can take follows, 

Mk) = (B4) 

Here gi{k) are real functions of k. 



Thus only the functions gi{k) remain undetermined in 
an arbitrary scattering process. Yet their low k expan- 
sions take a universal form, up to the coefficients of the 
expansion [28]. It is these coefficients which need to be 
calculated on a case by case basis (see Eq. (47) for the 
low k expansion of gi{k) = g-i{k)). 
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